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Abstract 

Bethe approximation is shown to violate Bravais lattices transla- 
tional invariance. A new scheme is then presented which goes over 
the one-site Weiss model yet preserving initial lattice symmetry. A 
mapping to a one-dimensional finite closed chain in an external field 
is obtained. Lattice topology determines the chain size. Using recent 
results in percolation, lattice connectivity between chains is argued to 
be g ( rf ~J)~ 2 where q is the coordination number and d is the space 
dimension. A new self-consistent mean-field equation of state is de- 
rived. Critical temperatures are thus calculated for a large variety 
of lattices and dimensions. Results are within a few percent of exact 
estimates. Moreover onset of phase transitions is found to occur in 
the range (d — l)q > 2. For the Ising hypercube it yields the Golden 
number limit d > . 
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1 Introduction 



The Weiss theory of ferromagnetism, also known as the mean-field theory, 
has been proposed almost a century ago [1]. Its major achievement is to 
provide a simple and universal model for collective phenomena and phase 
transitions. In most cases it yields phase diagrams which are qualitatively 
correct. However, quantitative results are very poor. In particular, all aspects 
of the critical behavior are grossly misrepresented. 

Thirty years latter, Bethe extended the Weiss one-spin treatment to a 
cluster of (<? + l) fluctuating spins where q is the coordination number [2]. For 
the one-dimensional Ising model, this scheme reproduces the exact results. At 
higher dimensions, though some non-universal critical parameters are better 
represented, the critical behavior stays classical. 

However it was recently pointed that within a Bethe scheme different 
site magnetizations are produced respectively in the (q + 1) cluster and its 
surrounding [3]. Therefore the Bethe scheme violates Bravais lattices trans- 
lational invariance. This symmetry violation which has been overlooked for 
decades, is indeed consistent with a well known result. Bethe equations of 
state correspond to the exact solution of a model on the deep interior of a 
Cayley tree (which by definition, has no translational symmetry and exhibits 
singularities at its surface [4]). 

In this paper we present a new self-consistent approach which embodies 
Bethe ideas, including fluctuations of a spin cluster, but preserves the lattice 
translational symmetry. 

The outline of this paper is as follows. A generalized lattice version of the 
usual mean field treatment of the Ising ferromagnet is presented in Section 
2. The underlying use of symmetry is thus explicited. Along the same 
lines, the Bethe scheme is revisited in Section 3. It shows the initial lattice 
symmetry is always broken. Section 4 is devoted to the presentation of a new 
self-consistent treatment which does preserves the symmetry. The lattice 
is divided into finite closed loops of fluctuating spins and mean-field spins 
respectively. The problem is thus mapped onto a one-dimensional closed 
fluctuating chain in an external field which allows an exact solution. The field 
is h = 5Jm where J is the coupling exchange, m is the mean-field average spin 
magnetisation and 5 is the chain connectivity. Using very recent results on 
percolation [5] it is argued that 5 = 9 J" where d is the space dimension. 
Critical temperatures are calculated in Section 5 for a large variety of lattices 
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at several dimensions. Discrepancies with available exact estimates are within 
few percent. In the last Section phase transitions are found to occur only in 
the range (d— l)q > 2 for d-dimensional Bravais lattices with coordination q. 
In the hypercubic ferromagnetic Ising case it gives the Golden number limit 
d> 1±&. 

2 The mean field scheme 

Mean-field theory is basically a one-site approach. However we present here a 
new generalized lattice scheme which enhances its interplay with symmetry. 
We illustrate it in the case of a ferromagnetic Ising system on a lattice with 
q nearest neighbors. The Hamiltonian is, 

H = -JY,S i S j . (1) 

where Si = ±1 is an Ising variable, and interactions J are restricted to 
nearest neighbors (nn). A mean- field treatment of Eq. (1) processes indeed 
in two steps. 

2.1 Breaking the lattice symmetry 

Two interpenetrated lattices A and B are created within a staggered symme- 
try. Thermal fluctuations are then ignored on the B-sublattice. Spins Sf are 
substituted to their common thermal average < Sf > denoted m =< Sf >. 
At this stage spins Sf are no longer equivalent to their g nn S^. The trans- 
lational lattice symmetry has been broken. Eq. (1) turns to, 

H MF = -JJ2Sf E <Sf>, (2) 

i&A jeB,j=l 

which is a decoupled one-site Hamiltonian. The partition function can now 
be calculated integrating over Sf = ±1. Magnetization at site i G A is then, 

< Sf >= tanh{Kgm} , (3) 

where K = f3J, (3 = ks is the Boltzmann constant and T the tempera- 
ture. 
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2.2 Restoring the symmetry 

At this stage symmetry has been artificially broken through the introduction 
of two interpenetrated sublattices A and B. The second step in the mean-field 
scheme is then to restore the initial lattice symmetry. It is done requiring the 
equality < Sf >= m which in turn produces the well known self-consistent 
mean-field equation, 

m = tanh{Kqm} . (4) 

Although fluctuations have been included only at one sublattice, the sys- 
tem is again homogeneous with the same magnetization m at every lattice 
site. The critical temperature T c is given by, 

KY = \. (5) 

Numerical estimates from Eq. (5) are rather poor. For instance, it is = 
0.25 for the square lattice (q = 4) instead of the exact value 0.4407 [6]. 
Moreover Kf = 0.50 at d = 1 (q = 2) whereas it should be T c = 0. 

3 The Bethe approximation revisited 

We now present an unusual view point about he Bethe scheme to enlighten 
its symmetry content. First a finite cluster is considered. Second the Bethe 
approximation is implemented starting from a lattice. On this basis it then 
is shown translational invariance is always lost in the process. 

3.1 A cluster in a surface field 

We consider a finite cluster of one spin Si coupled ferromagnetically to q 
nn spins Sj with an exchange J. The associated partition function yields 
< Si >=< Sj >= as expected for a finite spin system. 

An external field h is then applied but only on the surface shell of the 
cluster, i. e., only on the q spins Sj and not on the central spin Si. The 
Hamiltonian is, 

Hi^-JSij^Sj-hj^Sj. (6) 

3=1 3=1 
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A non-zero magnetization is induced in the cluster from h 7^ which is a 
symmetry breaking field. However since the field is not applied to the center 
spin, the cluster is no longer homogeneous, i. e., < Si >^< Sj > with both 
thermal averages being a function of h. At order one in a = f3h we have, 

< Si >= qt&nh(K)a , (7) 

and 

< Sj >= {1 + (q- l)tanh 2 (K)}a . (8) 

We can then ask which value of the surface field h, if any, makes the 
cluster homogeneous with a non-zero magnetization < Si >=< Sj >^ 0. It 
leads to the self-consistent equation in a, 

, coshfo; + K) . , 

a={q-l)\n — -) -(. 9 

cosh(a — K) 

Expanding the right-hand side of Eq. (9) around a = two solutions are 
obtained. First, 

h = k B T{3(q-l)(K-K c )} 1 * , (10) 
for K > K* where (tanhK c = 

Kf = \^- (id 

Second we have h = for K < . 

We thus have found that for each value of temperature K > , there ex- 
ists a well defined non-zero value of h which makes the cluster homogeneous, 
i. e., < Si >=< Sj >t^ 0. However for K < such a field does not exist, 
i. e., only h = is compatible with an homogeneous cluster. In this case 
< Si >=< Sj >= due to the cluster finiteness. Otherwise < Si >^< Sj >. 



3.2 The Bethe implementation 

The Bethe model starts from a lattice in which a (q + 1) cluster similar to 
the one treated above, is singled out. However the difference from above 
comes from the definition of the field h. Above, h is an external magnetic 
field. Here h is argued to represent an effective field which accounts for 
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the coupling between the cluster surface and the remainder of the lattice. 
Calculations are nevertheless exactly identical to those done in Subsection 
(3.1) leading to Eqs. (7) and (8). 

It is then usually stated that common magnetization in the cluster, i. e., 
m =< Si >=< Sj > is the magnetization in the lattice bulk. Eq. (7) gives, 

q 

m = a , (12) 

q-l 

where K is taken at K c . 

For a good review of the classical presentation of the Bethe model, see 
the textbook of Pathria [7]. Major improvement over Weiss model is twofold. 
Values of critical temperatures are slightly better. It is = 0.35 instead 
of 0.25 for the square case. At d — 1 results are exact with in particular 
K? = 0. 

3.3 The Bethe violation of translational invariance 

Nevertheless it has been proven recently [3] that a Bethe implementation on 
a lattice breaks translational invariance. To enlight this statement we present 
here a generalized lattice view similar to the one developped above for the 
mean-field model . 

To implement a Bethe scheme on a lattice 3 distinct interpenetrated sub- 
lattices must be introduced artificially. We start from one fluctuating spin 
labelled i (Si) to which we refer as an A-spin. The set of its nn spins are 
labelled by index j (Sj) and called B-spins. Last, nn spins of j spins, except 
the % spin, are labelled by the index k (S l ?) and called C-spins. The A-spin 
plus these B and C shells constitute a supercell for the lattice. This supercell 
is then used to pave the space and reproduce the whole lattice. Eq. (1) 
writes, 

H = -J"£{Si t S}+ ± Sj?}} . (13) 

ieA jeB,j=i jeB,j=i fcec* 

The prime over last sum means the summation over k j-na spins does not 
includes the i-spin. The scheme is easily visualized in the square case. For 
lattices with closed loops nn like triangular, Eq. (13) has to be modified to 
account for interactions between nn S*. 
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3.3.1 Breaking the lattice symmetry 

It is worth noticing Eq. (13) is identical to Eq. (1). At this stage no 
approximation was done. Only a formal relabelling has been used. 

Now we perform a mean-field transformation of the Hamiltonian. More 
remote C-sublattice spins S^ 3 are substituted to their thermal average mc =< 
S % k >. In so doing the lattice symmetry has been artificially broken. How- 
ever it allows a calculation of the partition function. Eq. (13) reduces to 
H = J2ieA Hi with, 

H^-JS^Sl-hj^S), (14) 

3=1 3=1 

where h = JJ2l=\ < S % k > an< ^ ^ runs f° r a ^ 1111 °f respectively each 
j spin, spin % being excluded. Hence 

h=J(q-l)m c , (15) 

for the effective field. From Eq. (14) we calculate thermal averages rriA =< 
Si > and niB =< S* > as a function of m c . We end up with three unknown 
magnetizations which are respectively rriA, and mc- 



3.3.2 The symmetry breaking is irreversible 

At this stage, as in the mean-field scheme, the artificially broken symmetry 
has to be restored to recover the initial physical problem we started from. 
It is of importance to stress even a drastic approximation should not change 
the symmetry of the problem. It means here the 3 different magnetizations 
have to be identical. 

First the symmetry within the fluctuating cluster A-B can be restored 
putting uia = tub- This equality yields a self-consistent equation in mc 
identical to Eq. (9), 

1 cosh(if [(g - l)m c + 1]) 

m c = — In — — — -T- , 16 

K cosh(A[(g — l)m c — 1J) 

where a = K(q — \)m c has been used. 

Next step is to extend the cluster homogeneity to the C-sublattice. De- 
noting m =< Si >=< Sj > we have to fulfill the equality m c = m. However 
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Eq. (7) shows that m = qtanh(K)K(q — l)mc which in turn obviously 
demonstrates m ^ mc, as long as Eq. (16) holds. 

Alternatively we could have achieved vtla = mc- Then it is at the expense 
of getting another self-consitent equation in mc, 

m c = q (q - 1)K t&nh(K)m c , (17) 

which is nevertheless not compatible with Eq. (16). In this case m\ ^ tub- 
In both cases, translational invariance is lost. There exist one super- 
sublattice A and B (A and C) embedded in the C-sublattice (B-sublattice), 
both having different magnetizations. We thus have clearly demonstrated 
that the Bethe scheme produces always an irreversible lattice symmetry 
breaking. 

3.3.3 The Bethe scheme is forbidden by symmetry 

Overlooked for several decades, this symmetry breaking could indeed be 
found readily from the Bethe scheme cluster character. Cluster center has 
all its nn spins which are fluctuating while surface cluster spins have (q — 1) 
mean-field nn spins and one fluctuating spin. Simultaneously mean-field 
spins have all their nn which are fluctuating spins, making their environ- 
nement identical to the cluster center. On this basis the Bethe requiremnt 
< Si >=< Sj > is not compatible with the equality m c =< Sj > which 
should also hold to ensure translational invariance. The Bethe topology is 
therefore forbidden by symmetry [3]. It is not the case for one-site Weiss 
theory. 

Last but not least, it is worth noticing it is this very symmetry problem 
which makes the Bethe model exact on the Caley tree lattice. This lattice 
does not exhibit translational invariance by construction. However the Cay- 
ley tree (also called Bethe lattice in the tree's deep interior) is pathological 
with singularities at its surface [4]. Spins on the surface have a non-negligible 
contribution to the partition function in the thermodynamic limit. This is 
not the case for a regular lattice. 

4 A preserving symmetry scheme 

Once the Bethe model is found forbidden by symmetry, the underlying ques- 
tion is to find out if a self-consistent treatment which includes more than one 
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site is indeed possible. We show below that such a model does exist. 

4.1 A 2-spins model 

Ogushi [8] considered a two nn fluctuating spin cluster. For one pair (1, 2) 
of fluctuating spins the equality < Si >=< S 2 > holds by construction. It 
allows then to fulfill < Si >= m as required by lattice symmetry. Critical 
temperatures are given by, 

(q-l)K°{l+tanh(K°)} = l . (18) 

Improvement over the Weiss model is only marginal with, for instance, K° = 
0.26 for the square lattice. Moreover a transition is still wrongly predicted 
at d — 1 . 

4.2 A TV-spin model 

To go over a two-spin cluster while preserving initial lattice symmetry, re- 
quires to start with a topology of only two spin species A and B. Therefore 
fluctuating spins Si make one species A and mean-field spins characterize 
the other species B. A two-sublattice coverage has to be built up under the 
constraint of identical A-spins and B-spins with respect to their respective 
topology (nn- surroundings) . 

4.2.1 Breaking the symmetry 

Above Bethe approximation analysis emphazises the role of the cluster center 
in the breaking of symmetry. Therefore to avoid such a symmetry breaking 
no fluctuating center should exist within the fluctuating cluster. Such a 
constraint can be achieved via compact closed linear loops within the lattice 
topology. For instance 4-spin squares and 3-spin triangles for respectively 
square and triangular lattices. 

Each one of these plaquettes is then set respectively as A-species and B- 
species with a staggered-like covering pattern. A-plaquettes (B-plaquettes) 
have thus all their nn plaquettes as B-plaquettes (A-plaquettes) . For a given 
plaquette, each spin has two nn spins of the same species within the plaquette 
itself and (q — 2) nn spins of the other species belonging to nn plaquettes. 
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At this stage the problem has been turned into to a series of fluctuating 
one-dimensional finite closed chains in an external field h. The number of 
spins N is determined by the lattice topology. It is N = 4, N = 3, N = 3 and 
N = 6 for respectively square, triangular, Kagome and Honeycomb lattices. 
Interactions with nn mean-field spin plaquettes produce the field h. We have 
h = 5Jm where 5 accounts for connectivety to B-sublattices. 

Above transformation allows then an exact solution since the problem of 
a nn Ising ferromagnet on a closed chain in an external field can be solved 
exactly [7]. In particular, the chain site magnetization is, 



< Si >= Pe^K ^-^^JJ lh , (19) 



;i - tanh(K) JV 
[1 + tanh(AT) 

at order one in h. Here i G A-plaquettes. 



4.2.2 Restoring the lattice symmetry 

Then putting < Si >= m restores initial lattice symmetry. Here symmetry 
restoring is possible since only two sublattices were involved. It was not the 
case for the three sublattice Bethe scheme. The self-consistent equation of 
state is, 

m = ^ex P 2 <^g) }m + ..., (20) 

at order one in m using h = 5Jm. To solve Eq. (20) we need to determine 
the value of S. 



4.2.3 Rescaling the connectivity 

At this stage it is worth evoking a recent work on percolation thresholds [5]. 
It was found that relevant connectivity variables for site and bond dilution 
are repectively (d — l)(q — 1) and I^zlKizll. In other words, starting from a 
given site, it shows that for site percolation, the number of possible directions 
(q — 1) has to be multiplied by (d — 1). For bond percolation this effective 
number of sites has to be divided by dimension d. Using these variables, a 
universal power law form was found to yield all percolation thresholds for all 
lattices at all dimensions [5]. 

This percolation finding suggests to consider here a rescaled connectivity 
between closed loops instead of 5 = q—2. Using above counting, we first start 
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with q instead of (q— 1) since now dealing with pair exchange interactions and 
not percolation. Then we renormalize q by (d — 1) giving q(d — 1). However 
the 2 neighboring sites which are treated exactly within the closed loop have 
to be substracted from the effective number of sites which gives q(d — 1) — 2. 
Moreover, interactions being related to bonds, we divide this number by d 
as for bond percolation which results in, 

s= ^-])-\ (21) 

(Jj 



5 Results 

We can now check the validity of our simple symmetry preserving model with 
respect to critical temperatures. From Eq. (21) we get, 

5K C exp 2K C { ^ - } - 1 • (22) 

The trivial connectivity counting 6 = q — 2 already improves Weiss model. 
For instance = 0.29 in the square case and T c — at d — 1. We now 
proceed using Eq. (21) for connectivty. 



5.1 d = 2 

For the square case (q = 4, N — 4), S — 1 which gives = 0.4399. Exact 
result is K e c = 0.4407. In the case of triangular lattice (q — 6, N — 3), 
= 0.2919 with 5 = 2 while the exact estimate is K e c = 0.2746. For 
Kagome (q — 4, = 3) 5 = 1 yielding = 0.4649 for an exact estimate of 
K c e = 0.4666. And = 0.6160 for the honeycomb lattice (q = 3, N = 6) 
where 5 = | for an exact estimate of K e c = 0.6585 (see Table I). 



5.2 d = 3 

Going to d — 3 imposes to restrict the plaquette size to = 4 since a one- 
dimensional loop cannot embody a three-dimensional topology. However 
there exits a d— dependence through Eq. (21). We get 6 = y, 5 = y and 
5 = y for respectively cubic, fee and bec lattices. Corresponding critical 
temperatures are given by = 0.2012, 0.1568, 0.1096 respectively for 
exact estimates of 0.2217, 0.1575, and 0.1021 (see Table I). 



11 



5.3 d>3 



Critical temperature estimates [9, 10] are available for the hypercube at d — 
5, 6, 7. These are K e c = 0.1139, 0.0923, 0.0777 respectively. In comparaison, 
Eq. (22) gives 0.1380, 0.0869, 0.0737. We can also predict values for other 
lattices like for instance fee ones (see Table I). 

5.4 d — > oo 

To get the d — > oo asymptotic limit of our model we take both q — > oo and 
J — > under the constraint gJ = est. From Eq. (21) connectivity limit is 
5 — > g(l — |) which gives always, 

5 -> g (23) 

at leading order. Indeed g diverges always quicker than ^.even for /cc- lattices 
where q = 2d(d— 1). In turn Eq. (22) becomes, 

K G C = - q , (24) 

which is exactly Eq. (5). Mean-field result is thus recovered as expected in 
the d — > oo limit. 



5.5 A^-dependence 

To evaluate the sensibility on the loop size, it is fruitful to expand Eq. (22) 
in powers of K. It gives 

K?(l + 2K G C + ... + ^^)(1 - {K?) N + ...)(1 " {K?) N + -..) = \ ■ (25) 
Since > 3, a simple analytic expression is obtained only at order one, 

K G C = \ . (26) 

At two dimensions Eq. (26) gives K G = 1, |, 1, 2 for respectively the square, 
triangular, Kagome and for honeycomb lattices. These results are rather poor 
and shows the importance of the finite value of which embodies part of 
the lattice topology. 
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5.6 Weiss model with 5 connectivity 



It is interesting to distinguish in our model the linear closed loop contribution 
from the use of the rescaled connectivity 5. It is readily done taking back 
Weiss Equation with g( - d ~ 1 - > instead of q. Here we don't have to substract the 
two sites of the linear loop as in Eq. (21). We thus get, 

K ? - th) < 27 > 

instead of Eq. (5). Critical temperatures at d = 2 are, K 9 C = 0.5000 for the 
square (q = 4, K e c = 0.4407), K% = 0.3333 for triangular lattice (q = 6, K e c = 
0.2746), K§ = 0.5000 for Kagome (q = 4, K e c = 0.4666) and K% = 0.6667 for 
the honeycomb lattice (q — 3, K e c = 0.6585) (see Table I). 

At d = 3 we get, K 9 C = 0.2500, 0.1875, 0.1250 respectively for exact 
estimates of 0.2217, 0.1575, and 0.1021 (see Table I). 



6 Conclusion 

We have presented a very simple self-consistent model which both preserves 
initial lattice symmetry and goes beyond the one-site Weiss approach. It 
yields values for critical temperatures within a few percent of exact results. 
Besides a rescaled lattice connectivity, the finite length of the loops is also 
taken into account. This new scheme represents a substantial improvement 
over existing mean-field cluster approximations. However associated critical 
exponents stay identical to those of the Weiss model. 

Last but not least, we can determine from our model a lower critical 
dimension for phase transitions. It comes from the condition h = for which 
we have a one-dimensional finite system. Such a system has no long range 
order at T ^ 0. Phase transitions are thus obtained only in the range h ^ 
which gives, 

q{d-l)>2. (28) 
For Ising hypercubes (q = 2d) it yields the Golden number limit, 

*>^, (29) 
which both excludes the d — 1 case and contains d — 2 as it should be. 
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Dimension 


Lattice 


Q 


5 






K% 


d = 2 


Square 


4 


1 


0.4407 


0.4399 


0.5000 




Honeycomb 


3 


2 


0.6585 


0.6160 


0.6667 




m • i 

triangular 


6 


2 


0.2746 


0.2837 


0.3333 


d = 2 


Kagome* 


4 


1 


0.4666 


0.4649 


0.5000 


d = 3 


Diamond 


4 


2 


0.3698 


0.2857 


0.3750 




sc 


6 


10 

3 


0.2216 


0.2012 


0.2500 




bcc 


8 


14 
3 


0.1574 


0.1568 


0.1875 




fee 


12 


22 

3 


0.1021 


0.1096 


0.1250 


d = 4 


sc 


8 


22 

4 


0.1497 


0.1380 


0.1667 




fee 


24 


23 
2 




0.0749 


0.0555 


d = 5 


sc 


10 


38 
5 


0.1139 


0.1064 


0.1250 




fee 


40 


158 
5 




0.0298 


0.0312 


d = Q 


sc 


12 


29 

3 


0.0923 


0.0869 


0.1000 


d = 7 


sc 


14 


82 

7 


0.0777 


0.0737 


0.0833 



Table 1: and K g c from this work compared to "exact estimates" K e c taken 
from [9, 10]. 
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